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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

A TRANSONIC INVESTIGATION BY THE FREE-FALL METHOD OF AN 
AIRPLANE CONFIGURATION HAVING SVfflPTBACK 
WING AND TAIL SURFACES 
By Stanley Faber and John M. Eggleston 

SUMMARY 


As part of a general investigation to obtain the characteristics 
of airplanes in the transonic range, the Langley Laboratory, by using 
the free-fall technique, has tested an aiiY>lane configtjration having 
45 ° sveptback wing and tail surfaces . The model was equipped with an 
all -movable horizontal tail that was stepped cyclically during the 
drop to vary the angle of attack. In this manner, the lift, drag, and 
longitudinal static and dynamic stability and control characteristics 
were obtained at various operating lift coefficients up to 0.26 over a 
Mach number range of O .5 to 1.21. The Reynolds number range was from 
2.5 X 10^ at release to 16 .O x 10^ at impact. Comparisons were made 
with test results of aerodynamlcally similar models of other trainsonic 
facilities in order that correlation between these facilities can be 
obtained . 

The resTilts of free-fall tests showed that only small trim changes 
in angle of attack and trim lift coefficient were experienced at high 
subsonic speeds . The lift-curve slope increased with Increasing Mach 
number up to a Mach number of 0.97 thereafter decreased. 

At Mach numbers near 1.2, the drag results indicate that some 
leading-edge suction was realized but that it was only a moderate per- 
centage of the amount predicted by theory. 

The ability of the horizontal tail to change the lift coefficient 
increased with increases in Mach number up to Mach number O .9 and there- 
after decreased. At the highest Mach nvmiber of the test, the tail was 
only half as effective in changing the trim lift coefficient as it was 
at low subsonic speeds . The model experienced a large increase in sta- 
bility at a Mach number of 0.9; the aerodynamic -center position moved 
from 40 percent of the mean aerodynamic chord rapidly rearward until at 
a Mach number of l.l? itie aerodynamic center was at 70 percent of the 
mean aerodynamic chord. 
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The longitudinal damping derivatives of the model indicated an 
increase in damping with increasing Mach number in the subsonic region 
and a more pronounced decrease in damping with increasing Mach nimiber 
supersonically . 

The stability parameters as calculated from the wind-tunnel data 
for pirrposes of comparison indicated the same general trends as the 
free-fall results . The most pronounced difference was in the Mach num- 
ber at which the large increase in static stability begins; the wind- 
tunnel results predicted this increase at a Mach number of approximately 
0.1 earlier than the free-fall results. Estimations of aeroelastic 
effects indicate that such effects were not a factor in these comparisons. 


INTRODUCTION 


As part of a general investigation to obtain the characteristics 
of airplanes in the transonic range, the Langley Laboratory, by using 
the free-fall technique, has tested an airplane configuration having 
45 ° sweptback wing and tail surfaces . The all -movable horizontal tail 
of this model was stepped cyclically in small increments during the 
drop in order to vary the angle of attack. In this manner, the transonic 
lift and drag, and the longitudinal static and dynamic stability and 
control characteristics were obtained at various operating lift coef- 
ficients (up to a value of 0.26) and over a Mach nimiber range of O .5 

to 1.21. The Reynolds number range was from 2.5 x 10^ at release to 

16 .0 X 10 ^ at impact . 

Directly comparable tests of aerodynamic ally similar models are 
being made by several other transonic facilities of the Langley 
Laboratory in order that correlation between transonic testing techniques 
can be obtained. 

The results of the free-fall test are presented herein as the vari- 
ation with Mach number of the lift coefficient, drag coefficient, lift- 
curve slope, pitching-moment slope, and damping-in-pitch derivatives. 

These results are compared with the parameters estimated from the tran- 
sonic wind-tunnel tests of the same wing -fuselage combination as reported 
in references 1 to 4 and with the zero-lift drag for the wing-fuselage 
combination as obtained from previous free-fall tests (ref . 5) • Frequency- 
response curves are also presented as obtained from the response of the 
model in angle of attack and normal acceleration to the step elevator 
deflections . 
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SYMBOLS 

exponential damping factor, as in e®*^, negative for stable 
oscillation, sec"^ 

mean aerodynamic chord, ft 

chord force coefficient, rW/qS 

lift coefficient, cos a - sin a 

incremental change in lift coefficient 
drag coefficient, Cjj sin a + cos oc 
pitching -moment coefficient, U/qSc 
normal-force coefficient, nW/qS 

variation of lift coefficient with angle of attack, per deg 
variation of lift coefficient with tail incidence, per deg 

variation of pitching -moment coefficient with angle of attack, 
per deg 

variation of pitching -moment coefficient with tail incidence, 
per deg 

variation of pitching -moment coefficient with pitching velocity 
( nondimens ional) ^ ^ 

variation of pitching -moment coefficient with rate of change 
of angle of attack ( nondimens ional)^ ^ 

~c d 

differential operator, 

" V dt 

frequency, cycles /sec 
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acceleration of gravity, ft/sec^ 

moment of inertia about lateral aocis, slug-ft^ 

tail incidence, deg 

incremental change in tail incidence, deg 
nondimens ional radius of gyration, l/mc^ 
tail length, ft 

distance from angle-of-attack vane to center of gravity, ft 
distance from normal accelerometer to center of gravity, ft 
mass, slugs 

Mach nvmiber; or pitching moment, ft-lb 
normal acceleration in g vmits 
static pressure, Ib/sq ft 

dynamic pressure, ^ pM^, Ib/sq ft 

longitudinal retardation in g units 

wing area, sq ft 

tail area, sq ft 

LaPlace transform operator 

time, sec 

velocity, ft /sec 

angle of attack, deg 

incremental change in angle of attack, deg 
rate of change of angle of attack, radians /sec 
ratio of specific heats, lA for air 


CONFIDENTIAL 


NACA RM L53D10 


CONFIDENTIAL 


5 


P 


de/da 


0 


0 


pitching velocity, radians/sec 
pitching acceleration, radians/sec^ 
variation of downvash angle with angle of attack 
airplane relative -density factor, m/pSc’ 
density, slugs/cu ft 


Subscripts ; 


f 


f\aselage 


t 


tail 


w 


wing 


The systems of a^es used in this report with positive directions and 
deflections are shown in figure 1. In plotting the frequency-response 
data, the sign of the tail incidence has been reversed to agree with 
dynamic -analys i s convent ions . 


The test model was a midwing airplane configuration as shown in the 
photograph in figiire 2. A drawing of the model giving the pertinent 
dimensions is presented as figinre 5 . 

The fuselage had a fineness ratio of 10, obtained by cutting off 
the rear one-sixth of a fineness -ratio -12 body. Table I gives the 
coordinates of the body. Reference 6 contains resiilts of drop tests 
of the body alone . The wing of the model had ^5° of sweepback measured 
at the quarter-chord line, an aspect ratio of A, a taper ratio of 0 . 6 , 
and an NACA 65 AOO 6 airfoil section parallel to the plane of symmetry. 

The wing had no dihedral or incidence and was located on the body such 
that the quarter chord of the mean aerodynamic chord was at 60 percent 
of the body length, the approximate model center of gravity. The all- 
movable horizontal tall was scaled down from the wing, the linear ratio 
being O.it-48, and was mounted in a slot in the vertical fin 26.4 percent 
of the wing semispan above the wing -chord plane . The pivot point of the 
tall was at a chordwise location of approximately 29 percent of the 
horizontal-tail mean aerodynamic chord. The slot in the vertical tail 
was sealed by cover plates attached to the horizontal tail. In order 
to keep the cover plates from blowing out, small shrouds were placed 
over the leading edge of the cover plates and were attached to the fin. 


MODEL 
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The junctiore of the tail surfaces incorporated a large rectangular 
fairing used primarily to provide the horizontal tail with adequate 
rigidity. The fairing also provided a flat sinrface for more effective 
sealing action of the cover plates. The fin had an NACA 65AOO9 airfoil 
section parallel to the body center line, an aspect ratio of 1.5^ a 
taper ratio of 0.5 and had the trailing edge swept back 45° • A photo- 
graph of the tail is shown in f ig\ire 4 . 

In order to obtain high strength and stiffness, the wing was made 
with a core of heat-treated steel. This core was covered with bismuth- 
tin alloy to form the actual airfoil shape . The vertical and horizontal 
tails were machined from solid aluminvun. 

The alinement of the model was very carefully controlled during 
construction in order to keep the rate of roll low. The necessity for 
keeping the rate of roll low stems from the desire of keeping the effects 
of rolling on the model stability small, as discussed in reference 7. 

Ground measurements showed the wing to have zero twist and an incidence 
angle of less than 1 minute . The vertical tail had a slight twist of 
about 4 minutes over its span, which was estimated to produce a rate of 
roll of one revolution per minute . 

The model weighed IO5O pounds with the center of gravity at 25.64 per- 
cent of the wing mean aerodynamic chord. The loading of the model was 
heavily concentrated near the model center of gravity in order to reduce 
any gyroscopic and centrif\;igal moments that would be produced by rolling. 

The moment of inertia of the model about the lateral axis throxigh the 
center of gravity was 50.I slug-ft^ and the wing loading was 114.5 Ib/sq ft. 

The incidence of the horizontal tail was changed diaring the drop in 
small increments by a stepping cam which was designed to give 2 seconds 
of fixed incidence at each position. The measured incidence angles were 
1.15°, 0.52°, -0.12°, -0.60°, -1.12°, -1.58°, -2.07°, and -2.87°. A 
time delay was employed to prevent the tail mechanism from operating 
during the initial I5 seconds so that the model could increase in speed 
at essentially zero lift and so that the model would not be at high lift 
coefficients near a Mach number of 0.9* Preliminary investigations 
showed the longitudinal stability to be marginal at this Mach number . 


MEASUREMENTS 


The test was performed by utilizing the free-fall method in which 
the flight path of the freely falling test body is obtained by radar and 
phototheodolite equipment and the other desired quantities are measured 
at the model by means of the NACA radio -telemetering system. References 5 
and 6 contain more complete details of this technique . An atmospheric 
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survey of static pressure^ temperature, and wind direction and velocity 
was made after the drop. The quantities measured and recorded by the 
telemetering system were: longitudinal retardation, normal acceleration, 

and transverse acceleration, all by accelerometers located 22 inches 
ahead of the center of gravity; the angle of attack, by a delta-shaped 
vane 7 inches ahead of the model nose (see fig. 5) i ihe tail deflection, 
at the root of the horizontal tail; the ambient pressiare, inside the 
base of the conical body of the angle -of -at tack vane; the impact pressure, 
by a probe located 5 inches back from the body nose and 1 inch from the 
body sinrface (see fig. 5); and. the rate of roll, by a rate gyro mounted 
in the model nose . 

The quantities which were used to determine the Mach number during 
the drop are presented as time histories in figure 6. The Mach number 
was obtained by first differentiating the flight path with respect to 
time to obtain the velocity of the model relative to the ground. This 
velocity was combined with the wind velocity to obtain the true airspeed. 
The Mach number was then determined from the true airspeed and the 
atmospheric temperature . The Mach number could be obtained from the 
radar -phototheodolite data for all but the last portion of the drop; a 
ground haze obscured the model from the optical trackers after k6 seconds 
from release . An auxiliary measurement of the Mach number was obtained 
from the telemetered impact pressure and the survey static pressure. 

These data were used to extrapolate the Mach mmber variation to impact 
as shown in figure 6. The Mach number variation thus obtained is believed 
to be accurate to within ±0.01. 


REDUCTION OF DATA 


The acceleration and retardation data were reduced to normal- and 
chord-force coefficients by using the model weight and wing area and the 
appropriate Mach number and static pressure . Lift and drag coefficients 

and Cj) were obtained from the normal- and chord-force coefficients and 

the measiared angle of attack. The variation of lift coefficient with angle 
of attack was obtained from incremental changes in the trim value of the 
parameters and a at each stabilizer step. These values were corrected 

to the case of constant tail incidence through use of the formula 

Sm Ai-j~ 
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where Ait/A>: is the change in tail incidence divided by the change in 

the trim value of the angle of attack. The tail lift-curve slope was 
estimated from tunnel tests of the wing inasmuch as the wing and tail 
of the model were geometrically similar. The effects of aeroelasticity 
on the lift-curve slope of the tail of the free-fall model were found 
to be about the same as those for the wing used in the wind-timnel tests . 
Therefore, no corrections for aeroelasticity were made. 

Corrected values of the lift-curve slope were also obtained from 
the alternate formula 


The derivation of this eqxiation assimies the lift-curve slope of the 
wing and of the tail to be the same . Somewhat different values of the 
lift-curve slopes would be expected because of the effects of aeroelasticity 
as the wing and tail were made of different materials; however, as will 
be shown subsequently, this effect of aeroelasticity is small and consider- 
ation of such effects in the formula does not significantly change the 
final result. Downwash values for this correction and for all calculations 
in this report were obtained from wind-tunnel data of reference k for the 
tail location of the free-fall model. The two methods used to correct the 
model lift-curve slope gave essentially the same results. 

The rigid-wing lift-ciorve slopes were obtained by correcting the 
measured by the methods of reference 8 as applied to a solid wing. 

Reference 9 was also used to obtain the rigid wing Ct^ and gave results 

■KZw 


in good agreement with those of reference 8. 

The zero-lift drag coefficient at supersonic speeds was obtained 
from the formula 


This formula applies when no leading-edge suction is present. The zero-, 
lift drag coefficient was also estimated at M = 1.2 by using the test 
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data obtained in the Mach range 1.17 to 1.21. In this Mach number 
range was fomd to be 0.2l6. 

The variation of the pitching-moment coefficient with angle of 
attack Cir^ was calculated by use of the eqxiation 


-(2nf)% 

n J 


The effects on the frequency of the damping in pitch and the degree of 
freedom involving motion along the vertical axis were neglected because 
these effects were found negligible for this model. 

The damping coefficients + Cmq were calculated from the 
formula 






qSc® 


where V is the true airspeed, and a is the damping exponent. The 
derivation of the equation for the damping coefficients was based on 
the assumptions that the model had two degrees of freedom. 

The estimated maximum uncertainties of the basic coefficients 
obtained from the telemetered measurements are presented in the fol- 
lowing table for two Mach n-umbers. The vincertainties are based on an 
instrument error of ±1 percent of the instrument range . 


M 


^D 

c„ 

“a 

0.8 

± 0 .029 

±0 .0082 

±0 .0010 

1.2 

±.010 

±.0027 

±.0013 
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The maximum uncertainties of other pertinent quantities that do not 
vary vith Mach number are as follows; 


Angle of attack, deg ±0.15 

Tail incidence (measured at root) , deg ±0 .02 

Rate of roll, radians/sec ±0.02 

Transverse acceleration, g lonits i0.02 


The measiired response of the model to the step elevator deflections 
were analyzed by using the Prony method to obtain the transfer functions 
in angle of attack and normal acceleration. The Prony method is described 
in reference 10 and, in brief, this method of analysis assumes that the 
control input and transient response may be expressed by a series of 
exponentials. The coefficients of the resulting analytical expressions 
are solved by a least-squares method, and the transfer ftinction is then 
established by taking the LaPlace transform of these analytical expressions . 
Transfer functions were obtained from the response in angle of attack as 
recorded at a point 5.5 feet ( Iq) ahead of the center of gravity and from 
the response in normal acceleration as measured at a point 2.165 feet ( Zp) 
ahead of the center of gravity. With these transfer functions, the fol- 
lowing relationships were solved simultaneously to establish the transfer 
functions in angle of attack and normal acceleration at the center of 
gravity. 
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These equations are based on the system of axes shown in figure 1. 
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RESULTS AND DISCUSSION 


General . - A full-size reproduction of a section of the telemeter 
record is shown in figure 7* Th^ high-frequency oscillation (approxi- 
mately ^7 cps) in the normal acceleration following the step is the 
first bending mode of the horizontal tail. Since the tail is swept back 
this bending produces oscillating tail loads which are reflected in the 
normal acceleration. The recorded amplitude of the oscillation is 
exaggerated somewhat by the accelerometer being located ahead of the 
center of gravity where it is sensitive to pitching accelerations. For 
the calculations based on these data this oscillation was faired. 

Complete time histories of the basic measurements obtained from 
the test are presented in figure 8. The figure shows that the time 
delay in the tail-programming device operated to produce the desired 
15-second time delay before the first step in control deflection occurred. 

The Mach number obtained at the time of this first step was 0.7. The 
programmer completed two cycles during the drop. At 4 seconds before 
impact^ the activating mechanism became jammed during the return stroke 
that finishes the cycle . This jamming was probably due to the hinge 
moments of the tail exceeding the output of the driving mechanism. The 
rate of roll can be seen to be small throughout the drop so that its 
effect on the stability of the model can be neglected. The transverse 
acceleration also was small; the small residual oscillation in the trans- 
verse acceleration which persisted throughout the drop is evidence that 
the configirration has some snaking tendencies. If a side-force coef- 
ficient of 0.01 per degree is assimied, the snaking oscillation had a 
fairly constant double amplitude of one -fifth of a degree. 

Lift characteristics .- The value of lift coefficient existing during 
the drop is shown in figure 9 with the corresponding incidence angles of 
the tail indicated. Because these lift coefficients were measured by an 
accelerometer located ahead of the center of gravity, the peak values of 
lift coefficient are affected by pitching acceleration. The oscillations 
shown were faired over small ranges of Mach number to obtain the average 
lift coefficient for a given tail incidence and Mach number. It should 
be remembered that variations with Mach number of the parameters presented 
herein should be associated with the particular trim lift coefficient 
which existed at each test Mach number. Examination of figure 9 shows 
quite clearly the effect of Mach number on the trim lift coefficients. 

Up to a Mach number of 0.91^ Mach number has little effect on the trim Cl 
for a given tail incidence but, above this value, the trim Cl decreases 
with increasing Mach number \mtil a Mach number of O.98, at which time a 
sharp increase in trim Cl occiors . An increase in trim Cl with 
increasing Mach number is also indicated in the Mach number range 1.125 
to 1.15. The maximum lift coefficients reached were 0.26 in the Mach num- 
ber range of I.05 to I.08 and 0.255 in the Mach number range of I.165 to I.I85. 
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The lift-CTjrve slope of the complete model is shown in fig- 


ure 10 and is compared with the 


calculated from the wind-tunnel 


tests of references 1 to 4 and in the subsonic region with the theoretical 
C]^ from reference 11. Also shown are the lift-curve slopes of the 

equivalent rigid wings . The lift-curve slopes were taken at the value of 
the mean lift coefficient existing at each test Mach number. The lift- 
curve slopes for the wing-body alone from the wind-tunnel data were cor- 
rected to the complete airplane configuration through use of the expression 





The theoretical wing lift-curve slope was corrected to the complete 
airplane lift-curve slope by the same expression. The values of 

determined from the free-fall test increased from the subsonic value of 
0.07 to a maximiam of O.IO8 at approximately M = 0.97 then decreased 

to a value of O.O7 at M = 1.2. At subsonic speeds the variation of CIq, 

obtained from the wind-tunnel tests has the same general shape as from 
the free-fall test, but the values from the wind-tunnel tests are smaller 
by a fairly constant amount of 0 .OO5 ♦ In the high subsonic range , the 
tunnel CIq does not peak as sharply as was indicated by the free-fall 

tests . As the Mach number is increased above unity, the decrease in the 
tunnel Ci^ is not so great as that shown by the present test, the tunnel 

value eventually becoming greater than free-fall value at the highest test 
Mach numbers . The theoretical curve has the same values as the wind- 
tunnel results for Mach numbers below 0.8 but does not indicate the sharp 
increase with increasing Mach number that both the free-fall and the wind- 
tunnel resiilts show. 

In order to investigate the effects of aeroelasticity on the comparison 
of the free-fall and wind-tunnel results, the lift-curve slopes for the 
rigid wing were calculated by using the methods of reference 8. For the 
free-fall data, the wing stiffness was based on that of the steel core of 
the model wing. For the timnel data, reference 3 shows that, in the range 
of lift coefficients of the present tests, no measurable differences in 
lift coefficient were found between a steel and an aluminimi wing. The 
steel wing was used for calculating aeroelastic effects on the wind-tunnel 
data. Correcting for the effects of aeroelasticity of the free-fall model 
tended to increase the value of the rigid-wing lift-curve slope over that 
of the flexible wing. The magnitude of this effect increased throughout 
the drop primarily because of the associated increase in dynamic pressure . 
The increase in lift-curve slope was negligible at M = O.5 and. increased 
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to about l6 percent at M = 1.2. The wind-tunnel tests, with a smaller 
variation in dynamic pressure, show an increase In rigid-wing lift-curve 
slope over that for the flexible wing ranging from 2 percent at low Mach 
numbers to 8 percent near the speed of sound and 7 percent at the highest 
Mach numbers . The comparison of the two models on the basis of rigid- 
wing lift-curve slope gives the same general conclusions as were shown 
by the measured results . 

Drag characteristics .- The variation of the drag coefficient Cp 
with Mach number is shown in figure 11. Also shown in the figure is a 
fairing of the zero -lift drag coefficient of the model at supersonic 
speeds estimated from the test data by using a lift-drag variation 
inversely proportional to the lift-curve slope. The zero-lift drag 
coefficient was also obtained at M = 1.2 through use of the lift-drag 
polar established from the test data in the Mach number range 1.17 "to 
1.21 and also thro\ 3 gh use of the methods of reference 12. As a basis 
for comparison, wind-tunnel drag coefficients corresponding to the mean 
lift coefficients at each stabilizer step of the present test and free- 
fall zero -lift drag from reference 5 are shown, both sets of data having 
been adjusted for the differences in configuration. For the wind-tunnel 
data, this adjustment involved addition of the estimated drags of the 
horizontal and vertical tails, correction of the base pressure drag, 
and correction for the sting error. Adjustments to the zero -lift drag 
of reference 5 included addition of the drags of the vertical and hori- 
zontal tail and correction for the difference in afterbody shape. Because 
of the differences in instrumentation the drag coefficients of the present 
test are only one-third as accurate as those of reference 5 - 

Figure 11 shows that, in general, the drag coefficients of all the 
tests are in good agreement. In the high subsonic range, the wind- 
tunnel resiats indicate an earlier drag rise, but the differences are 
within the estimated uncertainties in Mach number. At supersonic speeds, 
the results of reference 5 indicate a gradual increase in zero -lift drag 
coefficient with increasing Mach number. The estimated zero -lift drag 
coefficient of the present test also shows an increase but of a much 
reduced magnitude. From figixre 9, it can be seen that in a number of 
instances the operating lift coefficient was nearly the same over 
two ranges of Mach number. Comparison of the corresponding drag coef- 
ficients in figure 11 indicates that, under lifting conditions, there is 
an increase in drag coefficient with increases in Mach number . The rate 
of increase thus obtained is compatible with that of reference 5 « 

The zero -lift drag coefficient at M = 1.2 obtained by extrapolation 
of the measured lift-drag polar lies between the value estimated through 
use of the relationship involving the inverse of the lift-curve slope (no 
leading-edge suction) and the values estimated through use of reference 12 
(full leading-edge suction) . The position of the value obtained from the 
lift-drag polar indicates that some leading-edge suction was realized but 
only a moderate percentage of the f\ill theoretical amount . 
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Control characteristics .- The ability of the tail to change the 
model trim lift coefficient, plotted as Z^Xi/Ait, is shown in figure 12. 
The parameter increases from a subsonic value of 0.092 to a maximimi 
of 0.126 at M = 0.9 and then decreases to a value of 0.045 8-t M = 1.2. 
The first part of the decrease, that just after M = O. 9 , was due to a 
gradiaal increase in the stability of the model due to the rearward 
shift of the aerodynamic center of the wing -body combination. Further 
decreases up to the maxlmijm Mach number of the test were due in part 
to this cause and in pairt to the decreases in the lift-cvirve slopes of 
the wing and tail as indicated by figure 10 . 

The variation of model lift coefficient with tail incidence is 
shown in figure 15 for a Mach nimiber of approximately 1.2. The vari- 
ation is shown to be fairly linear and shows that the trends in 

discussed above are not affected by lift coefficient in the range of 
lift coefficients of this test. This nearly linear variation of Z^C^/aI^ 

also validates the use of the parameter dit/da in the equation for 
determining the llft-ciirve slope. 

Static stability characteristics .- From the angle of attack and the 
normal acceleration records (see fig. 7 ) > "the frequency of the oscilla- 
tion following the step inputs were obtained. Figirre l4 shows these 
data converted to the variation in pitching -moment coefficient with angle 
of attack Cni^, about the model center of gravity. Also shown is Cn^ 

calculated from the wind-tunnel res\ilts, which was obtained by adding 
the experimental Cidq for the wing plus fuselage at the test Cl, cor- 
rected to the model center-of -gravity location, to the estimated 

contribution from the tail. The tail contribution was obtained from the 
expression 


Also presented in figure l4 are the values of obtained from both 

the free-fall and the wind-tunnel models corrected for the effects of 
aeroelasticity by the methods of reference 8. In the subsonic region, 
the Cn^ of the free-fall model had a slight variation around a mean 

value of -O.OI 6 up to M = O. 9 , where there was a stable break, 

reaching a maximum of -0.0595 at M = 1.13- The stability then decreases 
to a value of approximately -O.O 52 at M = 1.2. The values estimated 
from the wind-tunnel tests are in fair agreement with those of the present 
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test with the exception that the stable break In the pitching -moment 
curve occixrs at roughly O.IM earlier in the case of the wind-tunnel 
test, Slotted-throat data (ref. 2) were used to obtain the wind- 
tunnel value of however, the closed-throat data (ref. 1) would 

give essentially the same Mach number of the pitching-moment break. The 
data corrected to the rigid configuration show the same results as did 
the data of the flexible models, the value of for the rigid case 

being only slightly more stable than Cm^ for the flexible case by an 

amount roughly proportional to the dynamic pressvire. 

Figvire 15 shows the variation of the aerodynamic -center position 
with Mach n\mber of the complete configuration; again, a comparison is 
made with the wind-tunnel estimates. Also shown is the aerodynamic- 
center position of the rigid configuration. Subsonically, the aerodynamic- 
center position tends to move forward with increasing Mach number from 
48 percent mean aerodynamic chord to 4l percent mean aerodynamic chord 
at M = 0.9 and thereafter the configuration becomes increasingly more 
stable; the aerodynamic -center position reaching a value of approximately 
70 percent mean aerodynamic chord at M = 1.17- The wind-tvmnel data 
agreed with the test data at low Mach numbers and at M = 1.2. In the 
transonic range, the wind-tunnel data indicated a more rearward location 
of the aerodynamic center than did the free-fall data. 

The effects of aeroelasticity shown in flg\ire 15 for the case of 
the free-fall model resulted primarily from a greater reduction in the 
lift-curve slope of the tail than of the wing; this effect was due to 
the difference in material used for these components . This factor did 
not contribute to the effects of aeroelasticity calculated for the wind- 
tunnel results since the tail was assumed to be of the same material as 
the wing. The other factor which affects the shift in model aerodynamic- 
center position due to aeroelasticity is the shift in the aerodynamic 
center of the wing itself. This factor was estimated to be small for 
both models. The over-all effect on the free-fall model was that the 
rigid model would have a more rearward location of the aerodynamic center 
from that of the flexible model, the displacement being roughly proportional 
to the dynamic pressxire. 

Damping characteristics .- The damping of the model, shown in coef- 
ficient form is presented in figure 16 . The data have wide 

scatter due to the relatively small amplitude of the pitch oscillations. 
These small amplitudes were difficiilt to measure and also showed, in 
some instances, evidences of nonlinear damping characteristics. The 
general trend of the data indicated that the dan 5 )ing increased as the 
Mach number was increased to unity and thereafter decreased with further 
increases in Mach mmber. The damping in the transonic region was approxi- 
mately 50 percent greater than the damping at the low- and high-speed ends 
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of the test curve, M = O .5 and 1.2. The data shown in figure 16 for 
comparison were calculated from the wind-tunnel res\ilts hy assuming all 
the damping to be contributed by the horizontal tail. This dan^jing 
was obtained from the expression 


The estimates thus obtained are generally smaller than those for the 
free-fall data and any discrepancy may be due in part to the damping 
supplied by the wing -body combination. 

In order to evaluate the period and damping of the test airplane 
configuration with respect to the current handling -qualities specifica- 
tions (ref. 13 )^ the aerodynamic derivatives of the present test were 
used to calculate the dynamic characteristics of a full-size fighter 
airplane flying at ^,000 feet. Values for the example airplane of the 
quantities affecting the damping are t 5 q>ical of current design practices . 
These values are as follows: wing loading of 65 Ib/ft^, moment of inertia 

in pitch of 23,600 slug-ft^, and a wing span of 30 ft. As compared with 
the model constants, the wing loading was halved, the radius of gyration 
was increased by a factor of approximately J, and the size was increased 
by a factor of The calculations were made for the center -of -gravity 
pcjsition of the test model and for a center -of -gravity position that 
would give a minimum subsonic static margin of 5 percent of the mean 
aerodynamic chord. These resiilts are shown in figijre I 7 plotted as the 
cycles to damp to l/lO amplitude. The specification states that the 
oscillation shall damp to l/lO amplitude in one cycle . The data in fig- 
ure 17 show that for the rearward center -of -gravity location the speci- 
fication is met for all Mach numbers below unity, the cycles required 
being just less than one in the low subsonic range, 6 / 1 O of a cycle at 

a Mach number of 0.9^ and thereafter increasing to 2^ cycles at M = 1.2. 

For the model center -of -gravity location the curve was displaced such 
that only near M = O .9 was the specification approximated, the airplane 

requiring li cycles to damp to l/lO amplitude subsonically and 2^ cycles 

at M = 1.2. A fact of interest here is that present experience indicates 
the specification may be too severe . Many current airplanes do not meet 
this requirement at high altitudes . 

Frequency-response characteristics .- The foregoing discussion of 
period and damping characteristics of the test model describes its free 
motions. In order to fvirther describe the dynamic characteristics of 
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the model in response to elevator inputs^ the frequency response and 
associated transfer functions have been calculated from the test data. 

With the model considered to be a linear system at any one condition 
of altitude^ Mach number^ and so forth^ the frequency response in an^le 
of attack and normal acceleration at the center of gravity to horizontal- 
tail movement was obtained by using the Prony method of analysis (ref. 10) 
and assuming two degrees of freedom. These frequency responses at given 
values of Mach number and over a pertinent range of frequencies are pre- 
sented in figure l8 for angle of attack and in figure 19 for normal acceler 
ation. The transfer functions, which describe the inherent dynamic charac- 
teristics of the model in response to an arbitrary input are presented in 
table II together with the flight conditions at which they were obtained. 
The transfer functions may be derived from the longitudinal equations 
given in reference lU . Although the freely falling body was not kept at 
a constant Mach number or altitude, the records were analyzed over periods 
short enough so that it appeared justifiable to average changes in alti- 
tude and Mach number. Nevertheless, the records were long enough to give 
reasonable accirracy in analysis. The procedure used for obtaining the 
transfer functions was such as to compensate for the effects of initial 
conditions . 

Examination of figures l8 and 19 shows that the frequency at which 
the peaks in amplitude ratio occur increased with increasing Mach number. 
This is the trend expected for the usual subsonic airplane . The high 
peak values of amplitude ratio relative to their static values is indi- 
cative of the low damping of the model. At these low values of damping 
the natural frequency of the model may be considered to be nearly equal 
to the frequency at which these peaks in the amplitude ratio occur. 

Since, for the free-fall model, increasing Mach number was obtained 
at decreasing altitude, the effect of altitude cannot be completely 
separated from the effect of Mach number. It should be expected, however, 
that damping would increase with decreasing altitude. (This effect would 
be indicated by a reduction in peak an^litude ratio relative to the static 
value.) Figures l8 and 19 indicate that, with increasing Mach number and 
decreasing altitude, the damping in angle of attack and normal acceler- 
ation remain almost unchanged but with a slight tendency towards less 
damping. It would appear, therefore, that the reduction in damping 
associated with the increasing Mach nimiber as discussed previously tended 
to cancel out the improved damping obtained by operation at lower altitudes 

The trend of the amplitude -ratio curves for normal acceleration indi- 
cate an increase in magnitude over the entire frequency range as Mach 
number was increased and altitude decreased. This trend appears reasonable 
for the case where Mach number effects are relatively small because, imder 
conditions where the stability derivatives do not change appreciably, the 
magnitude of the curves increase as a function of the square of the air- 
speed when the altitude is held constant. Similarly, for the case of 
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constant altitude and with no Mach number effects on the stability 
derivatives, the amplitude ratios in angle of attack would be expected 
to originate at approximately the same value of static sensitivity and 
to reach peak values of approximately the same magnitude but would be 
shifted over the frequency range an amount proportional to airspeed. 

These trends are roughly borne out for the test model up to a Mach 
number of 0.9* Above this Mach number, the over-all reduction of the 
amplitudes of the frequency response in angle of attack primarily reflects 
ein increase in static stability of the model. The trend is actually 
greater than indicated since the larger Mach numbers were obtained at 
lower altitudes and the altitude effect should improve the frequency 
response . A resume of altitude effects on the transfer-function coef- 
ficients together with their effects on the frequency-response curves 
is given in reference I 5 . 

The response obtained at a Mach number of O .725 appeaxs to be some- 
what different from those at other subsonic Mach mmibers inasmuch as a 
higher degree of damping and generally lower amplitude ratios are indi- 
cated thro\:ighout the frequency range investigated. No irregularity was 
found in the flight records that wo\xld tend to dispute or disprove the 
data obtained at this Mach number. 


CONCLUSIONS 


As part of a general investigation to obtain the characteristics 
of airplanes in the transonic-speed range, a model of an airplane con- 
figuration having 45° sweptback wing and tail siirfaces has been tested 
by the free-fall technique over a Mach number range of O .5 to 1.21 in 
order to determine the lift, drag, and longitudinal stability and control 
characteristics. The conclusions from these tests are as follows: 

1. The configuration experienced only small trim cheinges in angles 
of attack and trim lift coefficient in the high subsonic range. 

2. The value of the lift-c\urve slope of the free-fall model increased 
appreciably with increases in Mach mmiber at subsonic speeds . The lift- 
curve slope peaked at a Mach number of 0.97 attained a value 50 per- 
cent greater than the value obtained at low subsonic speeds and at a 
Mach mjmber of 1.2. At transonic speeds, the increase in lift-curve slope 
with increasing Mach number was much greater than that predicted by the 
simplified lifting-surface theory. 

3 . The results of the drag tests indicate near a Mach number of 1.2 
that some leading edge suction was realized but that only a moderate per- 
centage of the amount predicted by theory. 
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4. The ability of the horizontal tail to change the lift coefficient 
increased with increasing Mach number up to a Mach number of 0.9 and 
thereafter decreased. At the highest Mach number of the test, the tail 
was only half as effective in changing the trim lift coefficient as it 
was at low subsonic speeds . This decrease at supersonic speeds was due 

to the reduction in lift-ciirve slope of the tail and to the large increase 
in static stability of the model. 

5 . The static stability of the free-fall model decreased slightly 
with increasing Mach number up to M = O. 9 , where the aerodynamic center 
was at i<-0 percent of the mean aerodynamic chord. As the Mach number was 
further increased, the static stability increased rapidly until at a 
Mach number of 1.17 Ibe aerodynamic center was at 7 O percent of the mean 
aerodynamic chord . 

6 . The longitudinal damping derivative of the model 

indicated an increase in damping with increasing Mach number in the sub- 
sonic region and a more pronounced decrease in damping with increasing 
Mach number supersonically. The maximimi damping was approximately 
50 percent greater than the high- and low-speed values . When the aero- 
dynamic derivatives of the model were used to calciilate the dynamic 
characteristics of a full-size fighter airplane with typical values of 
center-of-gravity position, moment of inertia in pitch, and wing load.ing, 
and flying at an altitude of y),000 feet, the damping response of the 
airplane met the handling -qualities specifications at subsonic Mach 
numbers but required a slightly greater number of cycles to damp to the 
required amplitude at supersonic Mach numbers . 

7 . The data calculated from wind-tunnel results for pxrrposes of 
comparison indicate the same general trends as the free-fall results , 
the most pronounced difference being in the Mach mrniber at which the 
large change in the static stability began. 

8. The frequency-response curves obtained for the test model show 
trends typical of swept-wing configurations operating at transonic 

speeds and, in general, bear out the trends of damping, static sensitivity, 
and static stability given in conclusions 5 to 7- 
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TABLE I 


COORDINATES OF THE FINENESS -RATIO -10 BODY 


in. 

Y, 

in. 

X, 

in. 

Y, 

in. 

0 

0 

48.00 

4.876 

. 6 o 

.277 

54.00 

4.971 

.90 

.558 

60.00 

5.000 

1.50 

.514 

66 .00 

4.955 

3.00 

.866 

72.00 

4.828 

6.00 

1.446 

78.00 

4.610 

9.00 

1.956 

84.00 

4.274 

12.00 

2.565 

90.00 

5.75^^ 

18.00 

3.112 

96.00 

3.051 

24.00 

5.708 

100 .00 

2.485 

30.00 

4.158 



36.00 

4.489 



42.00 

4.719 
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TABLE n 

TRANSFER FUNCTIONS 


g 

it 


-Ct 


2n 


D + 




8\X^y^ 


D2 + 



D + 

-Cmg^Ci^ - ^^LCIn^ 

1 

> 
1 


8(i2Ky2 


— - Y1 
iT" “ 


•2 

gc 


2^l 


+ 


5 . 

^q) 

D + 


■it 

+ 

8H%y2 

An%y^ 


2k/ci^ - % - Cn^ 

T) + 




AMKy2 

xJ t 

1 

<< 

ro 

1 ‘ 



Mach 

number 

Altitude, 

ft 

g/it, deg/deg 

n/it^ g/deg 

0.725 

32,000 

-4.8557s + 80.547 

-0.0599s^ - 0.02866s + 8.855 

s2 + 2.08s + 49.209 

s 2 + 2.08s + 49.209 

.775 

50,500 

-7.1520s + 90.672 

- 0 . 04675 ®^ " 0 . 22644 s + 18.565 

s2 + 1.25s + 62.40 

s^ + 1.25s + 62.40 

.840 

28,500 

-1.7597s + 115.854 

-0 .0572532 - 0.67041s + 28.048 

s2 + 1.825s + 79.69 

s2 + 1.823s + 79.69 

.890 

27,000 

-1.9125s + 195 '996 

- 0 . 04150 s 2 - 0.12528s + 46.055 

s 2 + 2 . 52 s + 99.99 

s2 + 2.52s + 99.99 

.950 

25,500 

- 0 .l 4 ls + 174.872 

-0.0636s2 + 0 . 8400 s + 60.920 

s^ + 2.55s + 161.92 

s2 + 2.55s + 161.92 

1.040 

22,000 

-0.990s + 587.08 

- 0 . 11546 s 2 - 1.3550s + 159.021 

8 ^ + 4 . 00 s + 571.0 

s 2 + 4 . 00 s + 571 -0 

1.170 

15,200 

5.1517s + 515.56 

-0.176732 + 0.5169s + 261.970 

s 2 + 4.30s + 555.5 

s 2 + 4.50s + 535.5 
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Figure 1.- Sketch of system of 
directions and 


axes used in this report with positive 
deflections indicated. 
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Figure 2.- Photograph of con^ilete model. 
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Tail Length 
Wing Area 
Tail Area 
Aspect Ratio 
Taper Ratio 


44.53 in 
9sgf 
LTZsqff 
4 

0.6 


Airfoil Secfion^A/ing $ tail 
NACA 65AC06 parallel 
to body cxis 

Airfoil Section^vertical toil 
NACA 65 AGO 9 porallel 
to body axis 


Figure 5«- Details and dimensions of complete model. Coordinates of body 
surface are given in table I. Dimensions are in inches or as noted. 
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Figure 4.- Photograph of horizontal and vertical tails. 
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Figure 5.- Photograph of angle-of-attack vane and total head pressiare pickup. 
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Figure 6 .- Q^aniities used 1:0 determine Mach number. 
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Figure 7.- Full-size reproduction of a section of the record film. 
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Figure 8.- 


Tlme histories of recorded data. 
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Figure 9 -- Variation of model lift coefficient with Mach number. 
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Figure 10.- Variation of the lift-curve slope with Mach number. 
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Figure 11.- Variation of drag coefficient with Mach number. 
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Figure 12.- Variation with Mach number of the ability of the horizontal 
tail to change the model lift coefficient. 
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Tail incidence j deq 


Figure 15.- Variation of the lift coefficient with tail- incidence angle 

at Mach numbers near 1.2. 
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Mach number , M 

Figure Variation of the static longitudinal stability derivative Cm 

with Mach number. Center of gravity located at 23.6k percent 5. 
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Figure I 5 .- Variation of the aerodynamic- center position of complete 

configuration with Mach number. 
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Free fall 




Mach number of the number of cycles required 

nf r-i ■to <iaiiip to l/io amplitude for an airplane 

of fighter size and weight flying at 30,000 feet. 
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Figure l8.- Frequency-response data for the response of the angle of attack 

to the step elevator deflection. 
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Figure 19 -- Frequency-response data for the response of the normal 
acceleration to the step elevator deflection. 
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